Surface versus bulk Mott transition in Ca z Lai_ x V03 
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The Mott insulator LaVC>3 is known to become metallic if La is replaced by small concentrations 
of Ca. Since surface quasi-particle spectra of various perovskites are more strongly correlated than 
their bulk spectra, the coexistence of a metallic bulk and an insulating surface layer might be feasible. 
To investigate this possibility the dynamical mean field theory is used to evaluate the quasi-particle 
spectra in the bulk and at the surface of Ca^Lai-^VOs for various Ca concentrations. 



Significant surface contributions to quasi-particle spec- 
tra of strongly correlated transition metal oxides have re- 
cently been identified in photoemission measurements on 
several perovskite materials: SrRuC>3 Sr^RuC^ 0. Q. 
Ca a: Sr 1 _ x V03 HHH, and Ca^Lax.^VOa 0. The latter 
series is particularly interesting since LaVC>3 is insulating 
and becomes metallic upon doping with Ca. By varying 
the photon energy Maiti et al. Q were able to show that 
surface correlation features are more pronounced than in 
the bulk, in accord with the observations on the other 
perovskites cited above. Thus, the coherent peak at the 
Fermi level is less intense and the satellite peak below the 
V t 2g band is stronger than observed in the correspond- 
ing bulk photoemission spectra. Even more remarkably, 
Maiti et al. identified a Ca concentration range where 
the bulk has become metallic while the surface still ap- 
pears insulating, before, at even higher Ca doping, it also 
becomes metallic. Although sample preparation and ex- 
perimental resolution make the data analysis difficult, 
these results address a very interesting issue, namely the 
possible coexistence of different bulk and surface phases 
in semi-infinite systems. In the case of magnetism, for 
example, it is well known that at finite temperatures the 
magnetic moment can be reduced in the surface layers of 
a ferromagnet, but that both have the same Curie tem- 
perature unless the surface is artificially decoupled from 
the bulk 0. 

The identification surface-induced enhanced correla- 
tions is clearly important in order not to interprete them 
erroneously as bulk correlations. For instance, earlier 
photoemission data on Ca^Sri^VOa |4j appeared to be 
in conflict with various low-frequency measurements [t|. 
While the latter data showed that this perovskite series 
is metallic for all Ca concentrations, the photoemission 
spectra seemed to suggest that SrVC-3 is strongly corre- 
lated and that CaV03 is on the verge to a metal-insulator 
transition. Recent photoemission measurements 0, [f| 
covering a wide range of photon energies demonstrated 
instead that the bulk spectra of both materials are consis- 
tent with metallic behavior, but that correlation features 
are much more pronounced at the surface. While SrV03 
is a cubic perovskite, the oxygen octahedra in CaV03 are 
orthorhombically distorted. This leads to a slight reduc- 
tion in tig band width and a broadening of the van Hove 
singularity. The density of states at Ep, however, is only 
weakly affected. That the bulk properties of SrVC>3 and 



CaVC>3 should indeed be very similar was recently also 
found by Nekrasov et al. [13 in bulk electronic structure 
calculations combining the the local density approxima- 
tion (LDA) with the m any - body dynamical mean field 
theory (DMFT) [H[H[l|. 

As shown by Liebsch 01 using surface DMFT cal- 
culations, the enhancement of surface correlations can 
at least in part be understood in terms of the reduced 
atomic coordination at the surface and the concomittant 
narrowing of the local density of states. Thus, the ratio 
of local Coulomb energy and band width U /W is effec- 
tively enhanced at the surface. Accordingly, the coherent 
peak in surface quasi-particle spectra is less intense and 
the lower Hubbard band is enhanced compared to the 
bulk. Other mechanisms, such as additional band nar- 
rowing due to reconstruction and distortion of oxygen 
octahedra, or less efficient screening of electron-electron 
interactions at surfaces, might further increase the im- 
portance of correlations in surface photoemission spectra. 
The recent work on Ca x Sri_ x V03 therefore established 
that photoemission is consistent with low-frequency bulk 
probes as long as surface effects are properly taken into 
account. 

The same conclusion applies to Sr2Ru04 for which pre- 
vious photoemission spectra seemed to contradict bulk 
de Haas- van Alphen measurements E| Il5j . Recent ex- 
perimental and theoretical work [3lll6| proved, however, 
that this discrepancy can be resolved by taking into ac- 
count the lattice reconstruction at the surface of Sr2Ru04 
which leads to significant changes in the photoemission 
spectra. 

In the present work we investigate surface correlations 
in the perovskite series Ca^Sri-^VOs at Ca doping lev- 
els where metallic bulk and insulating surface behavior 
might coexist. Since doping with Ca leads only to minor 
structural changes and nearly constant Hubbard U [l7j . 
the modifications of the electronic properties are almost 
entirely caused by the degree of band filling. Accor ding 
to theoretical predictions by Potthoff and Nolting [l8| . 
a phase with metallic bulk and insulating surface prop- 
erties cannot occur unless electronic coupling between 
surface and bulk is assumed to be weak. These results 
were derived for the semi-infinite half-filled simple cu- 
bic s band, and by adopting a linearlized version of the 
DMFT in which low- and high-frequency excitations are 
essentially decoupled. It is therefore not clear to what ex- 
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FIG. 1: Solid curves: isotropic bulk density of states of cubic 
(a) LaV0 3 (d 2 ) and (b) Cao.5Lao.5VO3 (d 1 ' 5 ). Dashed curves: 
local density of d xz , yz states in the first layer (Ef = 0). 



tent they are applicable to realistic multi-band systems 
at arbitrary band filling. Here we calculate surface quasi- 
particle spectra by using a realistic layer-dependent lo- 
cal density of states appropriate for Ca.rSr1_.cVO3, com- 
bined with multi-band DMFT Quantum Monte Carlo 
(QMC) calculations. 

Fig. 1 shows the bulk and surface density of states for 
LaVC>3 and Cao.5Lao.5VO3. The main difference is the 
band filling. While the d 2 occupancy of the V t^g bands 
makes LaV03 insulating 0], Cao.5Lao.5VO3 with its d 15 
occupancy is metallic. In the bulk, the insulator-metal 
transition occurs at about x ~ 0.2 |5j. Since our main 
interest here is to study the effect of band filling on this 
transition in the bulk and at the surface we adopt for sim- 
plicity the ideal cubic symmetry and neglect the distor- 
tion of oxygen octahedra. According to the LDA-DMFT 
calculations for CaV03 ^3] these structural changes have 
little influence on the ti g quasi-particle spectra. The bulk 
density of states is therefore similar to that of SrV03. 
It exhibits the characteristic asymmetric shape resulting 
from the pronounced planar nature of th e t<?.„ states which 
is found for various peroskite materials |2Qj . 

In a cubic structure, the tig bands are degenerate. The 
surface lifts this degeneracy since only the intra-planar 
d xy states exhibit strong dispersion within the surface 
layer while d xz _ yz bands arc modified due to the reduced 
atomic coordination (the z-axis specifies the surface nor- 
mal) . The layer dependent local density of states is cal- 
culated using a tight-binding formalism for semi-infinite 



systems |2l|. Details concerning the application of this 
method to ti g bands are given Ref. |22|. As shown in 
Fig. 1, the spectral weight is shifted from the low- and 
high-energy regions to intermediate energies close to Ep. 
Thus, although the full band width coincides with the 
one in the bulk, the local density of d XZtVZ states in the 
surface layer is effectively narrowed. The surface distri- 
butions shown in Fig. 1 (a) and (b) are qualitatively sim- 
ilar. The small difference arises from slightly different 
surface potentials which are employed in the electronic 
structure calculation for semi-infinite Ca^Lai-^VOs in 
order to ensure charge neutrality. In the deeper layers, 
the local density of d XZiVZ states approaches the bulk t2 g 
density rather quickly, the main modification consisting 
in an oscillatory deviation which decreases with increas- 
ing distance from the surface. In contrast to the d xz , yz 
densities, the local density of d xy states in the surface 
layer is almost identical to the isotropic bulk density. 
This is a consequence of the pronounced planar character 
of the t2 g states. 

The narrowing of the local density states in the surface 
layer has an important influence on the quasi-particle 
distributions which we evaluate by using the multiband 
QMC-DMFT. Essentially we are dealing with coupled 
impurity problems where layer-dependent self-energies 
must be calculated and iterated until self-consistency is 
achieved. Since the full solution of this problem for multi- 
band systems is computationally not yet feasible, we ne- 
glect the coupling between layer baths and treat bulk 
and surface layers independently. In previous single-band 
DMFT calculations for surfaces this coupling was found 
to be small Enforcing charge neutrality at the sur- 
face presumably also diminishes the importance of inter- 
layer coupling between the impurity baths. The key in- 
put in the DMFT calculation is therefore the orbital- and 
layer-dependent local density of states. As a result of 
the one-electron band narrowing we can expect a similar 
enhancement of surface correlation effects as discussed 
previously for SrV0 3 and SrRu0 3 H3- 

As mentioned above, LaV03 is insulating, but hole 
doping by replacement of La via Ca induces an insulator- 
metal transition at rather low Ca concentrations. This 
can be seen clearly in the different shapes of the 
imaginary-time Green's function for bulk Ca^Lai-^VOs, 
as shown in Fig. 2 for several values of x. Note that the 
value of G(t) at t (3/2 is representative of the weight 
of the quasi-particle peak at Ep, while G{(3) gives the 
occupancy per spin-band. In the case of LaV03 (x = 0) 
already for U — 4.5. ..5.0 eV there is little weight near 
r ~ (3/2, indicating the reduction of intensity near Ep 
as U approaches the critical value for the Mott transi- 
tion. Of course, since the DMFT-QMC calculations are 
carried out at finite temperatures, G{(3/2) cannot ap- 
proach zero, i.e., a finite spectral density at Ep remains 
even in the nominally insulating phase. According to 
the results shown at finite Ca concentrations, however, 
G{(3/2) diminishes much less rapidly, indicating the in- 
cipient metallicity for x > 0. Thus, in the range of 
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FIG. 2: Quasi-particle Green's function as function of imag- 
inary time for bulk t^g states of Ca. x h&i- x VOz derived from 
DMFT (/3 = 8). G{(3/2) is representative of the quasi-particle 
weight at E F . (a) x = 0; (b) x = 0.1; (c) x = 0.2; (d) x = 0.5. 
Plotted are G(r) for U = 3.0 . . . 5.0 eV [up to 6 eV in (d)] in 
steps of 0.5 eV (from below); J = 0.7 eV. 



Coulomb energies where LaVC>3 is insulating, i.e., for 
U ps 5.0 eV, Ca^Lai-^VOs for x > 0.2 is metallic. We 
conclude from these calculations that within our model 
the borderline doping level between insulating and metal- 
lic behavior is about x ps 0.1. 

According to these results we consider x = 0.1 as the 
most promising Ca concentration at which the bulk is 
already close to metallicity while the surface layer might 
still be insulating. Fig. 3 shows bulk quasi-particle spec- 
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FIG. 3: Quasi-particle density of bulk t2 g states of 
Cao.1Lao.9VO3 (d 1 ' 9 ) for (a) U = 4.0 eV, (b) 4.5 eV and 
(c) 5.0 eV derived from DMFT (solid curves); J = 0.7 eV. 
Dashed curves: bare density of states. 



tra for Cao.1Lao.9VO3 for various Coulomb energies in 
the critical region of the metal-insulator transition. The 
analogous surface spectra are shown in Fig. 4. All spectra 
exhibit the three-peak structure characteristic of strongly 
correlated systems: a coherent peak at Ep and upper 
and lower Hubbard bands which appear as satellites in 
photoemission or inverse photoemission data. As in the 
systems studied previously , as a consequence of the 
effective narrowing of the d xz ^ yz subbands, correlations 
are more pronounced at the surface than in the bulk: the 
Hubbard peaks are more intense and the coherent peak is 
a narrower. Since the QMC-DMFT calculations are per- 
formed at about 1450 K (/3 = 8) the coherent peak per- 
sists even in the insulating phase. Nevertheless, in view of 
the temperature dependence of this feature known from 
one-band systems |llll23| the spectra in shown in Figs. 3 
and 4 suggest that at the surface the likelihood for a gap 
to open at Ep at low temperatures is slightly greater than 
in the bulk. In the bulk there is clearly more spectral 
weight near Ep compared to pure LaV03 (x = 0) [22| . 
reflecting the beginning metallic bulk behavior at x = 0.1 
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FIG. 4: Quasi-particle density of surface d xy states (solid 
curves) and d X z,yz states (dashed curves) of Cao.1Lao.9VO3 



(d l u ) for U = 4.0, 4.5 and 5.0 eV derived from DMFT; 
0.7 eV. Dotted curves: bare densities of states. 
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Ca doping concentration. The surface spectra shown in 
Fig. 4, on the other hand, indicate stronger tendency to- 
wards gap formation, in particular for U — 5 eV. 

In spite of the more pronounced surface correlations 
discussed above, it is not clear whether extrapolation 
to low temperatures would yield different Mott transi- 
tions for the bulk and the surface layer. Even though 
the spectral weight at Ep is consistently smaller at the 
surface, both N s (Ep) and Nb(Ep) could still vanish at 
the same critical U . This would be analogous to the 



behavior found in semi- infinite Heisenberg models: the 
magnetization near the surface is smaller than in the 
bulk, but the critical temperature for the transition to- 
wards paramagnetism is the same 0. Only if surface 
and bulk are strongly decoupled their critical behavior 
might be different. A careful examination of this question 
would also require to account for the coupling between 
the layer dependent 'impurity' baths which we have ne- 
glected for computational reasons. Presumably this cou- 
pling whould diminish the differences between bulk and 
surface spectra. 

So far we have only considered the effect of band nar- 
rowing on the d XZj y Z states which hybridize preferentially 
in planes normal to the surface. Surface reconstruction or 
distortion of oxygen octahedra presumably would reduce 
hopping also between d xy states in the first layer, thereby 
leading to additional enhancement of surface correla- 
tions. Moreover, less efficient electronic screening pro- 
cesses near the surface could increase the local Coulomb 
energy. Both effects would make the surface electron- 
ically different from the bulk. It is doubtful, however, 
whether this decoupling would be strong enough to give 
rise to separate bulk/surface Mott transitions. 

In summary, we have calculated bulk and surface 
quasi-particle distributions for Ca^Lai-^VOa using the 
dynamical mean field theory. At integer occupancy (d 2 
for x — 0) this perovskite material is a Mott insulator. 
At small Ca concentrations, however, the bulk material 
becomes metallic. As a result of the reduced atomic co- 
ordination at the surface, two of the t2 g subbands exhibit 
appreciable band narrowing, giving rise to stronger corre- 
lation features in the quasi-particle spectra. In the range 
of critical local Coulomb interactions for the Ca concen- 
tration i = 0.1 we have found at the surface a tendency 
for the formation of an excitation gap at a lower value of 
U than in the bulk. Whether this implies a coexistence 
of metallic bulk and insulating surface behavior is, how- 
ever, difficult to estimate at present because of the finite 
temperature used in the QMC-DMFT calculations. Also, 
additional surface effects on the electronic structure such 
as reconstruction, the charge state of the surface layer, 
and weaker screening of the local Coulomb energy should 
be considered. 
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